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Abstract: In these lectures, I shall discuss small x physics and the consequences 
of the high gluon density which arises as x decreases. I argue that an understanding 
t— i ■ of this problem would lead to knowledge of the high energy asymptotics of hadronic 

processes. The high gluon density should allow a first principles computation of 
: these asymptotics from QCD. 

in ■ 

m : 

O I LECTURE I: LOTS OF PROBLEMS 

On ■ 

On ; 

^ A Introduction 

ft, ■ I think we all believe that QCD describes hadronic physics. It has been tested 

in a variety of environments. For high energy short distance phenomena, perturba- 
tive QCD computations successfully confront experiment. In lattice Monte-Carlo 
computations, one gets a successful semi-quantitative description of hadronic spec- 
^ tra, and perhaps in the not too distant future one will obtain precise quantitative 

agreement. 

At present, however, all analytic computations and all precise QCD tests are 
limited to the small class of problems which correspond to short distance physics. 
Here there is some characteristic energy transfer scale E, and one uses asymptotic 
freedom, 

a s (E) -> (1) 

as E — > oo 

One question which we might ask is whether there are any non-perturbative 
"simple phenomena" which arise from QCD which are worthy of further effort. 
The questions I would ask before I would become interested in understanding such 
phenomena are 
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• Is the phenomenon simple and pervasive? 

• Is it reasonably plausible that one can understand the phenomena from first 
principles, and compute how it would appear in nature? 

I will in this lecture try to explain a wide class of phenomena in QCD which 
are pervasive, and appear to follow simple patterns. I will then try to explain why 
I believe that these phenomena can be simply understood within QCD. 

B Total Cross Sections at Asymptotic Energy 

Computing total cross section as E — > oo is one of the great unsolved problems 
of QCD. Unlike for processes which are computed in perturbation theory, it is 
not required that any energy transfer become large as the total collision energy 
E — > oo. Computing a total cross section for hadronic scattering therefore appears 
to be intrinsically non-perturbative. In the 60's and early 70's, Regge theory was 
extensively developed in an attempt to understand the total cross section. The 
results of this analysis were to my mind inconclusive, and certainly can not be 
claimed to be a first principles understanding from QCD. 

The total cross section for pp collisions is shown in Fig. 1. Typically, it is 
assumed that the total cross section grows as ln 2 E as E — > oo. This is the so called 
Froisart bound which corresponds to the maximal growth allowed by unitarity of 
the S matrix. Is this correct? Is the coefficient of ln 2 E universal for all hadronic 
precesses? Why is the unitarity limit saturated? Can we understand the total cross 
section from first principles in QCD? Is it understandable in weakly coupled QCD, 
or is it an intrinsically non-perturbative phenomenon? 

C How Are Particle Produced in High Energy Collisions? 

In Fig. 2, I plot the multiplicity of produced particles in pp and in pp collisions. 
The last three open circles correspond to the pp collisions with the multiplicity at 
zero energy subtracted. The remaining open circles to pp. The x's are pp collisions 
without the multiplicity at zero energy subtracted. Notice that the open circles 
fall on roughly the same curve. The implication is that whatever is causing the 
increase in multiplicity in these collisions may be from the same mechanism. 

The obvious question is can we compute N(E), the total multiplicity of pro- 
duced particles as a function of energy? 

At this point it is useful to develop some mathematical tools. I will intro- 
duce some useful kinematic variables: light cone coordinates. Let the light cone 
longitudinal momenta be 

p ± = ^=(E±p z ) (2) 



Note that the invariant dot product 
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FIGURE 1. The total cross section for pp collisions 
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FIGURE 2. Multiplicity of produced particles in pp and pp collisions 



p-q = p t -q t - p + q - p q + (3) 

and that 

p + p- = \{E 2 - p 2 z ) = \{p 2 T + m 2 ) = X -m% (4) 

This equation defines the transverse mass itlt- (Please note that my metric is the 
negative of that conventionally used in particle physics. An unfortunate conse- 
quence of my education. Students, please feel free to convert everything to your 
favorite metric.) 

Consider a collision in the center of mass frame as shown in Fig. 3. The right 
moving particle has p\ ~ \[2 \ p z | and p\ ~ m\f \ p z |. For the colliding 
particles = m pro j ect n e , that is because the transverse momentum is zero, the 
transverse mass equals the particle mass For particle 2, we have p% = Pi and 
P2 = Pi- 

If we define the Feynman a; of a produced pion as 

x = pt/pi 

then < x < 1. The rapidity of a pion is defined to be 

V = \ ln( yPi/Pn) = ^ln{2p +2 /m 2 T ) (6) 

For pions, the transverse mass includes the transverse momentum of the pion. 

The pion rapidity is always in the range —ycu < V < Ucm where ycu = 
ln{p + 1 m pro j ect u e ) All the pions are produced in a distribution of rapidities within 
this range. 

These definitions are useful, among other reasons, because of their simple prop- 
erties under longitudinal Lorentz boosts: p ± — > K ±1 p ± where n is a constant. Under 
boosts, the rapidity just changes by a constant. (Students, please check this rela- 
tionship for momenta under boosts.) 

A typical distribution of pions is shown in Fig. 4. It is convenient in the center 
of mass frame to think of the positive rapidity pions as somehow related to the right 
moving particle and the negative rapidity particles as related to the left moving 
particles. We define x = p + /Pp rojectUe and x' = p~ /p wojectae and use x for positive 
rapidity pions and x' for negative rapidity pions. 

Of course more than just pions are produced in high energy collisions. The 
variables we just presented easily generalize to these particles. 

Several theoretical issues arise in multiparticle production. Can we compute 
dN/dyl or even dN/dy at y = 0? How does the average transverse momentum 
of produced particles < p T > behave with energy? What is the ratio of produced 
strange/nonstrange, and corresponding rations of charm, top, bottom etc at y = 
as the center of mass energy approaches infinity? 

Does multiparticle production as E — > oo at y = become 



(5) 



FIGURE 3. The collisions of two ultrarelativistic nuclei as seen in the center of mass frame 
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FIGURE 4. A typical pion rapidity distributions for hadronic collisions. 

• Simple? 

• Understandable? 

• Computable? 

D Deep Inelastic Scattering 

In Fig. 5, a cartoon of deep inelastic scattering is shown. Here an electron 
emits a virtual photon which scatters from a quark in a hadron. The momentum 
and energy transfer of the electron is measured, and the results of the break up are 
not. In these lectures, we cannot develop the theory of deep inelastic scattering. 
Suffice it to say, that this measurement is sufficient at large momenta transfer Q 2 
to measure the distributions of quarks in a hadron. 

To describe the quark distributions, it is convenient to work in a reference 
frame where the hadron has a large longitudinal momentum p\ adron - The cor- 
responding light cone momentum of the constituent is vtonstituent- We define 
x = P constituent / P hadron- (This x variable is equal to the Bjorken x variable, which 
can be defined in a frame independent way. In this frame independent definition, 
x = Q 2 /2p ■ Q where p is the momentum of the hadronic target and Q is the 
momentum of the virtual photon. Students, please check that this is true.) The 
cross section which one extracts in deep inelastic scattering can be related to the 
distributions of quarks inside a hadron, dN/dx. 

It is useful to think about the distributions as a function of rapidity. We define 
this for deep inelastic scattering as 

V = y hadron -ln(l/x) (7) 

and the invariant rapidity distribution as 
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FIGURE 5. A cartoon of deep inelastic scattering. 

dN/dy = xdN/dx (8) 

In Fig. 6, a typical dN/dy distribution for a constituent gluons of a hadron 
is shown. This plot is similar to the rapidity distribution of produced particles in 
deep inelastic scattering. The main difference is that we have only half of the plot, 
corresponding to the left moving hadron in a collision in the center of mass frame. 

We shall later argue that there is in fact a relationship between the structure 
functions as measured in deep inelastic scattering and the rapidity distributions for 
particle production. We will argue that the gluon distribution function is in fact 
proportional to the pion rapidity distribution. 

The small x problem is that in experiments at Hera, the rapidity distribution 
function for quarks grows as the rapidity difference between the quark and the 
hadron grows. This growth appears to be more rapid than simply | y pro j — y | or 
{Vproj — y) 2 , and various theoretical models based on the original considerations of 
Lipatov and colleagues suggest it may grow as an exponential in | y pro j — y |. [1] 
If the rapidity distribution grew at most as y 2 , then there would be no small x 
problem. We shall try to explain the reasons for this later in this lecture. 

In Fig. 7, the Zeus data for the gluon structure function is shown. [2] I have 
plotted the structure function for Q 2 = 7 GeV 2 and 20 GeV 2 . The structure 
function depends upon the resolution of the probe, that is Q 2 . Note the rise of xg(x) 
at small x, this is the small x problem. I have also plotted the total multiplicity 
of produced particles in pp and pp collisions in the open circles on the same plot. 
Here I have used that y = log(E cm /l GeV) for the pion production data. This is 
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FIGURE 6. A rapidity distribution for gluons in the hadron wavefunction. 




FIGURE 7. The Zeus data for the gluon structure functions. Error bars on the data are not 
shown, but are about 10 percent 



FIGURE 8. A picture of a hadron viewed head on. 

approximately the maximal value of rapidity difference between centrally produced 
pions and the projectile rapidity. The total multiplicity has been rescaled so that 
at small x, it matches the gluon structure functions. (Strictly speaking, we should 
have plotted the total multiplicity at y — 0, but this is hard to extract from the 
data. If the distribution is an exponential in rapidity, then up to a constant these 
would be proportional.). Observe that the qualitative similarity between the gluon 
structure function and the total multiplicity. 

Why is the small x rise in the gluon distribution a problem? Consider Fig. 
8, where we view hadron head on. [3] The constituents are the valence quarks, 
shown as solid circles, and the gluons and sea quarks shown as circles with dashed 
lines. As we add more and more constituents, the hadron becomes more and more 
crowded. If we were to try to measure these constituents with say an elementary 
photon probe, as we do in deep inelastic scattering, we might expect that the 
hadron would become so crowded that we could not ignore the shadowing effects 
of constituents as we make the measurement. (Shadowing means that some of the 
partons are obscured by virtue of having another parton in front of them. For hard 
spheres, for example, this would result in a decrease of the scattering cross section 
relative to what is expected from incoherent independent scattering.) 

In fact, in deep inelastic scattering, we are measuring the cross section for 



a virtual photon 7* and a hadron, a^* hadron- Making x smaller correspond to 
increasing the energy of the interaction (at fixed Q 2 ). An exponential growth in 
the rapidity corresponds to power law growth in 1/x, which in turn implies power 
law growth with energy. This growth, if it continues forever, violates unitarity. 
The Froissart bound will allow at most In 2 (1/x). (The Froissart bound is a limit 
on how rapidly a total cross section can rise. It follows from the unitarity of the 
scattering matrix.) 

We shall later argue that in fact the distribution functions at fixed Q 2 do in fact 
saturate and cease growing so rapidly at high energy. The total number of gluons 
however demands a resolution scale, and we will see that the natural intrinsic scale 
is growing at smaller values of x, so that effectively, the total number of gluons 
within this intrinsic scale is always increasing. The quantity 

1 rfJV 

= ^ * < 9 > 

defines this intrinsic scale. Here nR 2 is the cross section for hadronic scattering 
from the hadron. For a nucleus, this is well defined. For a hadron, this is less 
certain, but certainly if the wavelngths of probes are small compared to R, this 
should be well defined. If 

A 2 »A 2 QCD (10) 

as the Hera data suggests, then we are dealing with weakly coupled QCD since 
a s (A) « 1. 

Even though QCD may be weakly coupled at small x, that does not mean 
the physics is perturbative. There are many examples of nonperturbative physics 
at weak coupling. An example is instantons in electroweak theory, which lead to 
the violation of baryon number. Another example is the atomic physics of highly 
charged nuclei. The electron propagates in the background of a strong nuclear 
Coulomb field, but on the other hand, the theory is weakly coupled and there is 
a systematic weak coupling expansion which allows for systematic computation of 
the properties of high Z (Z is the charge of the nucleus) atoms. 

If the theory is local in rapidity, then the only parameter which can determine 
the physics at that rapidity is A 2 . (Locality in rapidity means that there are not 
long range correlations in the hadronic wavefunction as a function of rapidity. In 
pion production, it is known that except for overall global conserved quantities 
such as energy and total charge, such correlations are of short range.) Note that 
if only A 2 determines the physics, then in an approximately scale invariant theory 
such as QCD, a typical transverse momentum of a constituent will also be of order 
A 2 . If A 2 >> 1/R 2 , where R is the radius of the hadron, then the finite size of the 
hadron becomes irrelevant. Therefore at small enough x, all hadrons become the 
same. The physics should only be controlled by A 2 . 

There should therefore be some equivalence between nuclei and say protons. 
When their A 2 values are the same, their physics should be the same. We can take 
an empirical parameterization of the gluon structure functions as 
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FIGURE 9. The Caldwell plot for the derivative of F 2 
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(11) 



7ri? 2 dy x 5 

where 5 ~ .2 — .3. This suggests that there should be the following correspondences: 

• RHIC with nuclei ~ Hera with protons 

• LHC with nuclei ~ Hera with nuclei 

To get some rough idea of what the scales are which are important, consider 
the Caldwell plot of the Zeus data shown in Fig. 9. [2] The function dF^j 'dlnQ 2 , the 
derivative of the structure function Fi is a function of both x and Q 2 . The Caldwell 
plot takes a slice of this data in the Q 2 ,x plane. The important thing to observe 
is the qualitative change in the behavior of the function at Q 2 ~ 3 GeV 2 at an 
x ~ 10" 4 . Until we get to this value, the curve is adequately parameterized using 
DGLAP evolution equations and GRV parameterization of the gluon distribution 
functions. It is suggestive that the turnover has something to do with the physics of 



high gluon density. If so the typical scale associated with gluon momenta is rather 
large. 

Since the physics of high gluon density is weak coupling we have the hope that 
we might be able to do a first principle calculation of 

• the gluon distribution function 

• the quark and heavy quark distribution functions 

• the intrinsic px distributions quarks and gluons 

We can also suggest a simple escape from unitarity arguments which suggest 
that the gluon distribution function must not grow at arbitrarily small x. The 
point is that at smaller x, we have larger A and correspondingly larger p T . A 
typical parton added to the hadron has a size of order 1/pt- Therefore although 
we are increasing the number of gluons, we do it by adding in more gluons of smaller 
and smaller size. A probes of size resolution Ax > 1/pr at fixed Q will not see 
partons smaller than this resolution size. They therefore do not contribute to the 
fixed Q 2 cross section, and there is no contradiction with unitarity. 



E Heavy Ion Collisions 

In Fig. 10, the standard lightcone cartoon of heavy ion collisions is shown. [4] To 
understand the figure, imagine we have two Lorentz contracted nuclei approaching 
one another at the speed of light. Since they are well localized, they can be thought 
of as sitting at x ± = 0, that is along the light cone, for t < At x ± = 0, the nuclei 
collide. To analyze this problem for t > 0, it is convenient to introduce a time 
variable which is Lorentz covariant under longitudinal boosts 

r = Vt 2 - z 2 (12) 

and a space-time rapidity variable 

"=5 fa (irf) (13) 

For free streaming particles 

z = vt = ^t (14) 
we see that the space-time rapidity equals the momentum space rapidity 

V = V (15) 

If we have distributions of particles which are slowly varying in rapidity, it 
should be a good approximation to take the distributions to be rapidity invariant. 




FIGURE 10. A space-time picture of ultrarelativistic nuclear collisions. 



This should be valid at very high energies in the central region. By the correspon- 
dence between space-time and momentum space rapidity, it is plausible therefore 
to assume that distributions are independent of rj. Therefore distributions are the 
same on lines of constant r, which is as shown in Fig. 10. At z = 0, r = t, so that 
r is a longitudinally Lorentz invariant time variable. 

We expect that at very late times, we have a free streaming gas of hadrons. 
These are the hadrons which eventually arrive at our detector. At some earlier 
time, these particle decouple from a dense gas of strongly interacting hadrons. As 
we proceed earlier in time, at some time there is a transition between a gas of 
hadrons and a plasma of quarks and gluons. This may be through a first order 
phase transition where the system might exist in a mixed phase for some length of 
time, or perhaps there is a continuous change in the properties of the system 

At some earlier time, the quarks and gluons of the quark-gluon plasma are 
formed. This is at some time of the order of a Fermi, perhaps as small as .1 Fermi. 
As they form, the particles scatter from one another, and this can be described 
using the methods of transport theory. At some later time they have thermalized, 
and the system can be approximately described using the methods of perfect fluid 
hydrodynamics. 

In the time between that for which the quarks and gluons have been formed 
and r = 0, the particles are being formed. This is where the initial conditions are 
made. 

In various levels of sophistication, one can compute the properties of matter 
made in heavy ion collisions at times later than the formation time. The problems 
are understood in principle for r > Tf ormation if perhaps not in fact. Very little is 
known about the initial conditions. 

In principal, understanding the initial conditions should be the simplest part 
of the problem. At the initial time, the degrees of freedom are most energetic and 
therefore one has the best chance to understand them using weak coupling methods 
in QCD. 

There are two separate classes of problems one has to understand for the initial 
conditions. First the two nuclei which are colliding are in single quantum mechani- 
cal states. Therefore for some early time, the degrees of freedom must be quantum 
mechanical. This means that 

AzA Pz > 1 (16) 

Therefore classical transport theory cannot describe the particle down to r = 
since classical transport theory assumes we know a distribution function f(p,x,t), 
which is a simultaneous function of momenta and coordinates. This can also be 
understood as a consequence of entropy. An initial quantum state has zero entropy. 
Once one describes things by classical distribution functions, entropy has been 
produced. Where did it come from? 

Another problem which must be understood is classical charge coherence. At 
very early time, we have a tremendously large number of particles packed into a 



longitudinal size scale of less than a fermi. This is due to the Lorentz contraction of 
the nuclei. We know that the particles cannot interact incoherently. For example, if 
we measure the field due to two opposite charge at a distance scale r large compared 
to their separation, we know the field fall as 1/r 2 , not 1/r. On the other hand, in 
cascade theory, interactions are taken into account by cross sections which involve 
matrix elements squared. There is no room for classical charge coherence. 

There are a whole variety of problems one can address in heavy ion collisions 
such 

• What is the equation of state of strongly interacting matter? 

• Is there a first order QCD phase transition? 

These issues and others would take us beyond the scope of these lectures. The 
issues which I would like to address are related to the determination of the ini- 
tial conditions, a problem which can hopefully be addressed using weak coupling 
methods in QCD. 



F Universality 

There are two separate formulations of universality which are important in 
understanding small x physics. 

The first is a weak universality. This is the statement that physics should only 
depend upon the variable 

A = (17) 

As discussed above, this universality has immediate experimental consequences 
which can be directly tested. 

The second is a strong universality which is meant in a statistical mechanical 
sense. At first sight it appear a formal idea with little relation to experiment. If 
it is however true, its consequences are very powerful and far reaching. What we 
shall mean by strong universality is that the effective action which describes small x 
distribution function is critical and at a fixed point of some renormalization group. 
This means that the behavior of correlation functions is given by universal critical 
exponents, and these universal critical exponents depend only on general properties 
of the theory such as the symmetries and dimensionality. 

Since the correlation functions determine the physics, this statement says that 
the physics is not determined by the details of the interactions, only by very general 
properties of the underlying theory! 

We can see how a renormalization group arises. In Fig 11, the space-time 
distribution of gluons is shown. The coordinate space rapidity is used. The effective 
action which we shall develop in later lectures is valid only for gluons with rapidity 
less than 770- [5]- [8] Those at larger rapidity have been integrated out of the theory 
and appear only as recoilless sources of color charge for r) < 77 < y pro j- 
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FIGURE 11. The space-time distribution function for glue inside of a hadron. 

The way a renormalization group is generated is by integrating out the gluon 
degrees of freedom in the range rji < rj < i] , to generate a new effective action 
for rapidity r\ < 7]\. We can show that this results in new source strength of color 
charge, now in the range rji < rj < Ufrag, and in a modification of some of the 
coefficients of the effective action. 

At high Q 2 , the renormalization group analysis simplifies, and one can show that 
in various limits reduces to the BFKL or DGLAP analysis. The renormalization 
group equations at smaller q 2 become more complicated, and have yet to be written 
in explicit form and evaluated. This is in principle possible to do. [9] 

An essential ingredient in this analysis is the appearance of an action and of 
a classical gluon field. We should expect the appearance of a classical gluon field 
when the phase space density of gluons becomes high. Gluons are after all bosons, 
and when the phase space density is large they should be described classically. 
This provides the essential difference between older renormalization group analysis 
which was formulated in terms of a distribution function. High density and its 
complications due to coherence require the introduction of a field. The new renor- 
malization group analysis is phrased in terms of the effective action for the classical 
gluon field and jumps from an ordinary equation to a functional equation. 

The classical gluon fields which we shall find are the non-abelian generalization 
of the Lienard-Wiechart potentials of electrodynamics. If we use the coordinate 
space variable x~ and realize that the source for the gluons arise from much smaller 
x~ than that at which we make the measurements, since they arise from gluons of 
much higher longitudinal momenta which are more Lorentz contracted, then the 
sources can be imagined as arising from a 5(x~). The Lienard-Wiechart potentials 
also are proportional to 5(x~), and so exist only in the sheet. The fields are also 
transversely polarized 



K J- K JL z 



(18) 



G Why an Effective Action? 

The effective action formalism which will be advocated in the next lectures is 
very powerful. It is used to compute a gluon effective field. This field can be related 
to the wavefunction of the hadron. 

This field allows one to generalize it from the original description of a single 
hadron, to collisions of hadrons and also to diffractive processes. We shall see that 
the effective action formulation is incredibly powerful. 

The careful student might at this point be very worried: How has the problem 
been in any way simplified? We have just introduced a source, and to specify 
the field one has to specify the source. Moreover, such a specification is gauge 
dependent. What happens is amusing: We integrate over the all color orientations 
of the source. Gauge invariance is restored by the integration. The theory becomes 
specified by the local density of color charge squared. This is a gauge invariant 
quantity. It can be related to the gluon distribution functions, and is determined 
by the renormalization group equations. 



II AN INTRODUCTION TO LIGHT CONE PHYSICS 

This lecture will provide an introduction to light cone kinematics and quanti- 
zation of field theory on the light cone. We will eventually use light cone methods 
to quantize QCD, using the light cone gauge. 

Light cone coordinates are 



x 



± 



1 



V2 

and momenta 



(x°±x 3 ) (19) 



P ± = ^(P°±P 3 ) (20) 

The invariant dot product is 

p • x — p t • x t — p + x~~ — p~x + (21) 

where p t and x t are transverse coordinates. This implies that in this basis the 
metric is g + ~ = g~ + = —1, g tj = 5 tj where i,j refer to transverse coordinates. All 
other elements of the metric vanish. 

An advantage of light cone coordinates is that if we do a Lorentz boost along the 
longitudinal direction with Lorentz gamma factor 7 = cosh(y) then p ± — ► e ±y p ± 



FIGURE 12. The initial data problem on the light cone. 



If we let x + be a time variable, we see that the variable p~ is to be interpreted as 
an energy. Therefore when we have a field theory, the component of the momentum 
operator P~ will be interpreted as the Hamiltonian. The remaining variables are 
to be thought of as momenta and spatial coordinates. In Fig. 12, there is a plot 
of the z, t plane. The line x + = provides a surface where initial data might be 
specified. Time evolution is in the direction normal to this surface. 

We see that an elementary wave equation 

(p 2 + M 2 )(j) = (22) 

is particularly simple in light cone gauge. Since p 2 = p\ — 2p + p~ this equation is 
of the form 

is first order in time. In light cone coordinates, the dynamics looks similar to that 
of the Schrodinger equation. The initial data to be specified is only the value of 
the field on the initial surface. 

In the conventional treatment of the Klein-Gordon field, one must specify the 
field and its first derivative (the momentum) on the initial surface. In the light 
cone coordinate, the field is sufficient and the field momentum is redundant. This 
means that the field momentum will not commute with the field on the initial time 
surface! 

Lets us work all this out with the example of the Klein Gordon field. The 
action for this theory is 



S = -jd'x[\{d^f+ l -M^^ 



(24) 



The field momentum is 

Note that II is a derivative of <fi on the initial time surface. It is therefore not an 
independent variable, as would be the case in the standard canonical quantization 
of the scalar field. 

We postulate the equal time commutation relation 

[n(x t ,x-),<f>(y t ,y-)\ = -i8W(x-y) (26) 

Here we time is x + = y + = in both the the field and field momentum. We see 
therefore that 

d_[0(£),0($] = -^ 3 )(x-y) (27) 



or 

[4>{x)A{y)\ = -ie(x~ - y-)5 (2 \x - y) (28) 

Here e(v) is 1/2 for v > and —1/2 for v < 0. 

These commutation relations may be realized by the field 

= l>o(2^{ e "' IaW + e ~"' VW } (29) 

Using 

[a(p).a*(q)] = 2p + (2nf 5 {z \p - q) (30) 

the student can verify that the equal time commutation relations for the field are 
satisfied. 

The quantity l/p + in the expression for the field in terms of creation and 
annihilation operators is singular when p + = 0. When we use a principle value 
prescription, we reproduce the form of the commutation relations postulated above 
with the factor of e(x~ —y~)- Different prescriptions correspond to different choices 
for the inversion of One possible prescription is the Liebrandt-Mandelstam 
prescription l/p + — p~ / (p + p~ +ie). This prescription has some advantages relative 
to the principle value prescription in that it maintains causality at intermediate 
stages of computations and the principle value prescription does not. In the end, 



for physical quantities, the choice of prescription cannot result in different results. 
Of course, in some schemes the computations may become prohibitively difficult. 

The student should now check that with the field above, the light cone Hamil- 
tonian is 



as it must be. 

In a general interacting theory, the Hamiltonian will of course be more com- 
plicated. The representation for the fields in terms of creation and annihilation 
operators will be the same as above. Note that all particles created by a creation 
operator have positive P + . Therefore, since the vacuum has P + = 0, there can 
be no particle content to the vacuum. It is a trivial state. Of course this must be 
wrong since the physical vacuum must contain condensates such as the one respon- 
sible for chiral symmetry restoration. It can be shown that such non-perturbative 
condensates arise in the P + = modes of the theory. We have not been careful in 
treating such modes. For perturbation theory, presumably to all orders, the above 
treatment is sufficient for our purposes. 



In QCD we have a vector field A%. This can be decomposed into longitudinal 
and transverse parts as 



and the transverse as lying in the tow dimensional plane orthogonal to the beam z 
axis. Light cone gauge is 




(31) 



A Light Cone Gauge QCD 




) 



(32) 



(33) 



In this gauge, the equation of motion 



D U F^ = 



(34) 



is for the + component 



DiF i+ - D + F- + = 



(35) 



which allows one to compute A in terms of A 1 as 



(36) 



This equation says that we can express the longitudinal field entirely in terms of the 
transverse degrees of freedom which are specified by the transverse fields entirely 
and explicitly. These degrees of freedom correspond to the two polarization states 
of the gluons. 

We therefore have 

- 4 « w = £>. wh ( e "" x(p) + e " x ' (p) ) (37) 

where 

K(P), a{\q)} = 2p + 5 ab ^{2nf5^\p - q) (38) 
where the commutator is at equal light cone time x + . 



B Distribution Functions 

We would like to explore some hadronic properties using light cone field oper- 
ators. For example, suppose we have a hadron and ask what is the gluon content 
of that hadron. Then we would compute 



dNgi uon 
d 3 p 



=< h | a ] (p)a(p) | h > (39) 



The quark distribution for quarks of flavor i (for the sum of quarks and antiquarks) 
would be given in terms of creation and annihilation operators for quarks as 

^ =< h | {b\(pM P ) + 4(p)d z (p)} | h > (40) 

where b corresponds to quarks and d to antiquarks. The creation and annihilation 
operators for quarks and gluons can be related to the quark coordinate space field 
operators by techniques similar to those above. The interested student should read 
the notes of Venugopalan for details. [10] 

How would we begin computing such distribution functions? We will start with 
the example of a large nucleus, as this makes some issues conceptually simpler. We 
will then generalize to hadrons, where we shall see that the ideas presented here 
have a generalization. 

For a large nucleus, we assume that the gluon distribution which we shall try to 
compute has longitudinal momentum soft compared to that of the valence quarks. 
Valence quarks have longitudinal momentum of the order of that of the nucleus, 
so this requirement is only that x « 1 In fact, we will require that x « A^ 1 ^. 
This is the requirement that in the frame where the longitudinal momentum of 
the gluons is zero, the nucleus has a Lorentz contracted size much less than the 
wavelength associated with the gluons transverse momentum A ~ This is 



dN 
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Units 

FIGURE 13. The rapidity distribution for quarks and gluons inside of a hadron. 

the requirement that the gluon resolve the nucleus as a whole and is insensitive to 
the details of the nuclear structure (spatial distribution of valence quarks within 
the nucleus). 

In Fig. 13, a rapidity distribution for gluons and for valence quarks is shown. 
We are interested in the region where the overlap between the quark and gluon 
rapidity distribution is small. 

If the gluon phase space density is very large, the quantum gluons of the nucleus 
may be treated as classical fields. This may be true if 

A 2 = 1 dN 9 lue (ai) 
7iR 2 dy K ' 

satisfies A 2 >> A.q CD . Certainly if the typical pr of the gluons was of order Aqcd 
this would be true since the gluons would then be closely packed together. We shall 
see that this is true in fact for gluons with p T « A which for high gluons density 
can become very large. 

If the valence quarks have a longitudinal momentum much larger than the 
typical gluon momentum then their typical interactions with these gluons should 
be characterized by the soft momentum scale (otherwise the soft gluons would not 
remain soft). In an emission of a gluon, the emitted gluon has momentum very 
small compared to the valence quark longitudinal momentum. Its velocity therefore 
is barely changed by this emission. The quarks are therefore recoilless sources of 
color charge for the gluon classical field. 

The picture we have is therefore that of Fig. 14. The valence quarks are 
sources for the gluons and sit on a sheet of thickness infinitesimal compared to the 
typical wavelength associated with the gluon field. We shall further assume for 
simplicity that the transverse distribution of charge is uniform. (In fact, a better 
approximation is that the distribution of charge is slowly varying on the gluon 
wavelength scale of interest. This case can be computed directly from knowledge 
of the uniform transverse distribution case, The transverse size scale of variation 
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FIGURE 14. The nucleus as a source of charge sitting on a Lorentz contracted sheet. 



A x « 1 Fm 



FIGURE 15. A head on picture of a nucleus. 



is the nuclear radius, which is much larger than a fermi, so that this criteria is 
satisfied so long as A >> Aqcd, and the typical transverse momentum is of order 
A.) 

In Fig. 15, a head on picture of the nucleus is shown. A small square is 
indicated which shows the transverse resolution size scale that we shall employ to 
probe the nucleus. In order that we can use weak coupling methods, we require 
that the transverse size be Ax « lFm. 

In Fig. 16, the nucleus is shown along the beam axis. Where the tube generated 
by the square associated with Ax intersects a quark or gluons shown an x. ( We 
use space-time rapidity variables to spread out the nucleus. The precise definition 
of this variable is below. For now, just imagine that we have chosen a longitudinal 
spatial variable which spreads out the Lorentz contracted hadron.) The physical 
extent in longitudinal spatial coordinate is of course small. If we require that 
Ax << lFm and require that a tube intersects a quark or gluon from some hadron, 
then typically, that hadron will be far separated from any other hadron which 




FIGURE 16. The nucleus as seen longitudinally with the longitudinal axis blown up by a choice 
of variable such as rapidity. 

has one of its quarks or gluons within the tube. The quarks and gluons within 
the tube are therefore uncorrelated in color. If we further require that the tube 
be large enough so that many quarks and gluons are in each tube, a situation 
always possible if we have large enough energy or a large enough nucleus, then the 
color total color charge associated with these quarks and gluons will be in a high 
dimensional representation, and can be treated classically. To see this, recall that 
in a high dimensional representation, Q 2 » Q ~ [Q, Q] so that commutators may 
be ignored. 

For partons at the x values where we wish to compute the gluon distribution 
function, the higher rapidity sources sit on a sheet of infinitesimal thickness. The 
problem we must solve is to compute the typical correlation functions which give 
distribution functions for stochastic sources sitting on a sheet of infinite transverse 
extent traveling at the speed of light. Mathematically, this problem is similar to 
spin glass problems in condensed matter physics. Such problems typically have 
entirely nontrivial renormalization groups, and we will see that this is the case for 
our problem. 

To understand a little better how the sources are distributed, it is useful to go 
inside the sheet. This is accomplished by introducing a space-time rapidity vari- 
able. Let us assume that the rapidity of the projectile is y pro j and its longitudinal 
momentum is P + , we define 

V = y P roj - ln(P+ roj x-) (42) 

We define the momentum space rapidity as 

y = - HP+ oj /p + ) (43) 

Here x~ is the coordinate of a source and p + its momentum. The previous definition 
of momentum space rapidity for produced particles was 



y = -ln(p + /p ) = ln(p + /m t ) (44) 

where m t = \Jp 2 + m 2 . This last expression is valid for particles on mass shell, 
whereas the other definition of momentum space rapidity works for constituents of 
the hadron wavefunction, which are of course not on mass shell. These definitions 
are equal within about a unit of rapidity for typical values of transverse momenta. 
In the wavefunction, the uncertainty principle relation gives x~p + ~ 1, so that 
we see the space time rapidity is up to about a unit of rapidity the same as both 
momentum space values. We will therefore use these rapidities interchangeably in 
what follows. 

We can now implement our static sources of charge by assuming the theory is 
defined by an ensemble of such charges: 



[dp]expl--J^ dy'd 2 x t ^—^p 2 {y\x t )\ 



(45) 



This ensemble gives 

< p a (y,x T )p b (y',x' t ) >= p 2 (y)5 ab 5(y - y')5 {2 \x t - x' t ) (46) 

The parameter p 2 therefore has the interpretation of a charge squared per unit 
rapidity 

My) = -1 ^-9M (47 ) 

M2/J N 2 -lnR 2 dy { ' } 
The total charge squared is 

X(V)= r rOJ dy'p\y') (48) 

Since the sources are individual quarks and gluons, this can also be related directly 
to the total number of gluons and quarks contributing to the source as [11] 

*<^i(t + l^)/> GW (49) 

The factors above arise from computing the color charge squared of a singe quark 
or a single gluon. 

We may now construct the non-abelian Lienard-Wiechart potentials generated 
by this distribution of sources. We must solve the equation 

D,F^ = g 2 5» + p(y, x t ) (50) 



To solve this equation, we look for a solution of the form 



A ± = 



(51) 



and 

A 1 = -U( y) x t )V l U\y,x T ) (52) 

We could equally well solve this equation in a gauge where A + is nonzero and all 
other components vanish. The student should check that the gauge transformation 
induced by U gives 

A + = -U ] d + U (53) 

i 

In this non-lightcone gauge, the field equations simplify and we get 

^A + = g 2 p(y,x t ) (54) 

Solving for U gives 

i J dy'^p(y',x t )\ (55) 

We have therefore constructed an explicit expression for the light cone field A 1 in 
terms of the source p for arbitrary p\ The system is integrable. 

The structure of the fields strengths E and B which follow from this field 
strength is now easy to understand. The only nonzero longitudinal derivative is 
d + . The field strength where both ij are transverse vanishes since the filed 
looks like a pure gauge transformation in the two dimensional space. It is in fact a 
pure gauge everywhere but in the sheet where the charge sits. Therefore the only 
nonvanishing field strength is F l+ . This gives E _L B _L z, that is the fields are 
transversely polarized. This is the precise analog of the Lienard-Wiechart potentials 
of electrodynamics. 



Ill LECTURE 3: THE ACTION AND 
RENORMALIZATION GROUP 

The effective action for the theory we have described must be gauge invariant 
and properly describe the dynamics in the presence of external sources, up to an 
overall gauge transformation which is constant in x ± . (The lack of precise gauge 
invariance arises from the the desire to define the intrinsic transverse momentum 
of gluon distribution functions. Although when used in computations of gauge 
dependent quantities, the gauge dependence disappears, it is useful to not have the 
full gauge invariance in intermediate steps of computation.) The student should 
verify, that consistency with the Yang-Mills equations in the presence of an external 
source requires that 



(56) 



For a source of the type we have here, 



J£ = 5» a 5(x-)p a (x t ) 



(57) 



This is an approximation we make when we describe the source on scales much 
larger than the spatial extent of the source. To properly regularize the delta func- 
tion, we need to spread the source out in x~ as was done in the previous section. 
We find that the action is 



In this equation, the matrix T is in the adjoint representation of the gauge group. 
This is required for reality of the action. The student should minimize this action 
to get the Yang-Mills equations, identify the current, and show that the current is 
covariantly conserved. 

This action is gauge invariant under gauge transformations which are required 
to be periodic in the time x + . This is a consequence of the gauge invariance of 
the measure of integration over the sources p. This will be taken as a boundary 
condition upon the theory. In general if we had not integrated over sources, one 
could not define a gauge invariant source, as gauge rotations would change the 
definition of the source. Here because the source is integrated over in a gauge 
invariant way, the problem does not arise. 

In the most general gauge invariant theory which we can write down is generated 
from 



This is a generalization of the Gaussian ansatz described in the previous lecture. 
It allows for a slightly more complicated structure of stochastic variation of the 
sources. The Gaussian ansatz can be shown to be valid when the evaluating struc- 
ture functions at large transverse momenta. 



This theory is a slight modification of what we described in the first lecture. 
We have here assumed that our theory is an effective theory valid only in a limited 
range of rapidity much less than the rapidity of the source. In this restricted range, 
the structure of the source in rapidity cannot be important, and therefore we couple 
only to the total charge seen at the rapidity of interest. The local charge density as 
a function of rapidity is replaced by the total charge at rapidities greater than that 
at which we measure the field. The scale of fluctuation of the source is instead of 
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(59) 




(60) 



the local charge squared per unit are per unit rapidity fi 2 (y) becomes replaces by 
the scale of fluctuation of the total charge. The student should prove that 

fVproj „ 

X = / dy'fiyi) (61) 

To fully determine F in the above equation demands a full solution of the 
renormalization group equations of the theory. This has yet to be done. We shall 
confine our attention in most of the analysis below to the Gaussian ansatz. It is 
remarkable that within this simple ansatz for F, most of the general feature which 
a full treatment should generate, such as proper unitary behavior of deep inelastic 
scattering, arises in a natural way. 

Let us turn our attention for the moment to the gluon distribution function 

(2tt) 3 2 P + ^ =< a\p)a{p) > (62) 
Using the results of the previous lecture, you should prove that 

Here a is a color index and % is a transverse index associated with the gluon field. 
D is the gluon propagator in the external field 

D%(p, q )=<AZ(p)A» b (q)> (64) 

(In field theory language, D is the propagator in the external field including both 
connected and disconnected pieces.) 

In lowest order, the A in the expression for D is simply the external Lienard- 
Wiechart potential. Recall that 

A 1 = -iUV'U* (65) 

where the C/'s are the explicit functions of the source p computed in the previous 
lecture. 

The expectation value above may be computed using the Gaussian weight func- 
tion. We need the propagator 1/V^ to perform this calculation. The student should 
try to do this computation (it is equivalent to normal ordering exponentials), and 
if it is too hard, refer to the paper of Jalilian-Marian et. al, where it is worked out 
in some detail. [7] The result is 

dN r,2^- iPTZT 4(iv c 2 -i) 



d 2 pq 



J d 2 ze~^ V " zthocn) {l - (z^lcn) 2 ™^} (66) 



This results in the form for the gluon distribution function as shown in Fig. 17. 
At p t » A where A = a sy /x, the 
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FIGURE 17. The intrinsic transverse momentum dependence of the gluon distribution function. 
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A 2 /a s p 2 t 



(67) 



This is because this part of the distribution can be thought of as arising from 
bremstrahlung from the sources of the gluon field at high rapidities. As pt < A, 
the gluon distribution function goes to a slowly varying function of px, which in 
the Gaussian ansatz is 



dN 



1 



ttR 2 dyd?pt 



~ — ln(A/p t ) 



(68) 



In general, we expect saturation in this region, for arguments which will be pre- 
sented in the next few pages. This means we get a slowly varying function of pr, 
which by dimensional arguments, is therefore a slowly varying function of A 

The only place that any non-trivial rapidity dependence enters the problem is 
through A. This dependence can be found from the renormalization group analysis 
of the effective action. Therefore in the bremstrahlung region, the dependence on 
A 2 is linear, and in the saturation region it is very weak. 

Let us assume that when we go beyond the Gaussian ansatz, the general features 
of saturation presented above remain. We will show that this results in a reasonable 
physical picture for small x physics, and solve the unitarity puzzle outlined in the 
first lecture. First consider the gluon distribution function. The total number of 
gluons which can be measured at some scale size larger than a resolution size scale 
1/Q 2 is 



xG{x,Q 2 ) 



Q 2 d 2 p T dN 
o (27r) 2 dyd 2 p T 



(69) 



For Q 2 » A 2 , the integral is dominated by the bremstrahlung tail, and G ~ 
R 2 A 2 . (The gluon distribution is proportional to the charge per unit area times 
the area.) In our random walk scenario, the effective charge squared must be 
proportional to the length of the random walk, R, so that G ~ R 3 which for 
a nucleus is proportional to A 1 / 3 . This is what we expect at large Q 2 , except 
the reasoning is a little different than usual. The standard argument would have 
been that at large Q 2 , the degrees of freedom in a nucleus for example should act 
incoherently, and we should have the gluon distribution functions proportional to 
A. Here we have random fields generating precisely the same A dependence! 

When Q 2 < A 2 , the integral is dominated by the saturation region. Here 
G ~ R 2 Q 2 , and the gluons can be thought of as arising from the surface of the 
hadron. Again this is consistent with what is expected from phenomenology. The 
gluons are so soft that they cannot see the entire hadron, only its surface. 

The effect of saturation for unitarization of deep inelastic scattering can also 
be easily understood. Suppose we are at some fixed Q 2 » A. 2 (y). As y decreases 
corresponding to going to smaller x,, the gluon distribution function increases as 
A. 2 (y). When x becomes so small that we get into the saturation region, the linear 
A dependence is weakened (in the Gaussian ansatz it is logarithmic in A), and the 
structure function stops growing. We expect that the number of gluons of size 
smaller than our resolution scale have ceased to grow, and the cross section should 
become slowly varying. This is in spite of the fact that A continues to grow! 

The physics is again simple to understand: At small x we are indeed adding 
more and more gluons to the hadron wavefunction. These gluons are smaller as 
their inverse size scales as gluon density, 

1 r 2 = — 70 

1 ttR 2 dy v ; 

They do not contribute to a fixed Q 2 cross section when they become smaller than 
the resolution size scale. What saves unitarity is p t broadening. 

A Deep Inelastic Scattering 

In the previous discussion, we were concerned with the gluon distribution func- 
tion. In deep inelastic scattering, we measure the quark distribution functions. 
How are these related? 

In deep inelastic scattering we measure 

W muu = ^_ Im J d 4 xe i gx < p had | T{r{x)r{0) | p had > 
1 f 
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P»-^)[p v -<f^)—F,} (71) 



We must be able to compute the imaginary part of the current- current correla- 
tion function. This is simply vacuum polarization in the presence of the Lienard- 
Wiechart potentials. In fact, it is straightforward to compute the propagators in 
these background fields. [13] The reason why it is simple is because the background 
field is basically gauge transformations in two regions of spaces separated by a 
surface of discontinuity. 

This propagator may be then used to compute the vacuum polarization. The 
computation can be done so far as to relate explicitly the structure functions to cor- 
relation functions of exponentials involving p and expressions very similar to those 
for the gluon distribution function result. These results are currently being applied 
to the study of deep inelastic scattering to see whether the effects of saturation 
might be seen experimentally. 

B Renormalization Group and How It Works 

In the above computation of the gluon propagator, one used the classical effec- 
tive action to compute the Lienard-Wiechart potentials. These classical fields were 
then used to compute the propogator. What about the quantum corrections? 

The classical field contribution is shown in Fig 18. The x in the figure marks the 
position of the source. In Fig 19, the contribution arising by inserting the lowest 
order contribution to the connected piece of the gluon propogator is shown. Here 
there is a loop digram and the dashed line with the solid ellipse in it represents the 
gluon propogator to all orders in the strength of the external field. In Fig. 20, the 
piece where a quantum loop correction to one of the classical fields is shown. 

When the quantum corrections are computed, one gets a correction to the 
classical field of order a s ln(x cuto ff/x), where x is the value for gluon distribution 
function, and x cuto ff is the maximum x appropriate for the effective action. Al- 
though the coupling is small, the logarithm can become big if we go to x values far 
below the cutoff. The overall theory should remain valid however, the problem is 
that the classical approximation is no longer good. 

The way to fix this up is to use the renormalization group. Suppose we want 
a new effective theory below x new , and we require that x new « x cuto ff but that 
a s ln(x cu tof f I Xnew) << 1- This means that the quantum corrections are small in the 
range x new << x << x cu t ff, and can be handled perturbatively. We proceed by 
integrating out the degrees of freedom in this intermediate range of x, to generate 
a new effective theory in the region x « x new . 

In this process, the one can show that the only thing that changes in the effective 
action formalism is the weight function F[p] for fluctuations in the color source . 
Effectively, we are changing what we call source for the gluon field, and trading it 
for what we call the gluon field. This renormalization group procedure is that of 
Wilson-Kadanoff. One can derive the form of the renormalization group equations 
and they can be written explicitly in the low gluon density region. One can prove 
that here the function F is Gaussian, and that in appropriate kinematic limits, 
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FIGURE 18. The classical field contribution to the gluon distribution function. 




FIGURE 20. The quantum propogator correction to the gluon distribution function. 
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FIGURE 21. The ln(l/x) and ln(Q 2 ) plane and various regions of high and low parton density. 

one reproduces the BFKL equation, the DGLAP equation and their non-linear 
generalizations to first order in the non-linearities. Work is currently in progress 
to get the explicit form of these equations in the high density region, and to solve 
them. 

The various regions of high and low density are shown in the plot of Fig. 21. If 
we are in the low density region, we evolve in x by the BFKL equation and in Q 2 
by the DGLAP equation. Hopefully a full solution to the problem will lead to an 
effective action which is at a fixed point of the renormalization group. In this case, 
at very small x, the form for F will simplify, and all correlation functions will have 
universal critical exponents. If this is true, the dynamics of high energy scattering 
for all hadrons becomes the same, and presumably has a simple structure. 

C Other Directions 

In Raju Venugopalan's talk, you heard about using these techniques to describe 
high energy hadron-hadron collisions. [14]- [16] I have nothing to add to his talk, 
except that this shows that there is a direct relationship between hadron-hadron 
collisions and deep inelastic scattering. 

A third class of phenomena is diffraction. In Fig. 22, a rapidity distribution 
of particle for a diffractive event is shown. There are two clusters of produced 
particles. In deep inelastic scattering these could be say the clump associated with 
a quark-anti-quark jet produced by the virtual photon, and a clump associated with 
the fragmentation of the target. If we restrict our attention to the class of diffractive 
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FIGURE 22. The rapidity distribution of particles produced in a diffractive event. 

phenomenon where the target does not fragment, that is scatters elastically, one 
can prove the following: [17]- [19] 

• Deep inelastic scattering is given by computing the amplitude for vacuum 
polarization for the electromagnetic current in the presence of the non-abelian 
classical field, squaring and then averaging over color. 

• Diffractive events are given by computing the amplitude in the presence of the 
non-abelian color field, averaging over color, and then squaring. 

We have argued the first case earlier in the lecture. To understand how the 
second case works, note that only the diagrams of Fig. 23 survive averaging over 
color before squaring the amplitude. In this figure a virtual photon produces a 
pair of quark jets. These jets exchange gluons (dash dot line) with the source 
which comprise the hadron. For sources are far away in rapidity from the jet, and 
therefore the longitudinal momentum transfer is very small. (The sources are not 
kicked by the jet. They do not change their velocity.) Moreover the color averaging 
and translational invariance guarantees that the transverse momentum given to the 
quarks after two scatterings vanishes. 

So if we imagine a hadron and its Fock space constituents, they are unchanged 
by the scatterings shown in Fig. 23. Therefore the initial hadron state projects 
on to the hadron after scattering, in the amplitude, with weight very close to one, 
that is the hadron has scattered elastically (although many of its constituents have 
scattered.) 

The amplitude shown in Fig. 23 therefore describes the rapidity distribution of 
particles shown in Fig. 22 only if the hadron does not fragment, that is, maintains 
its identity as a hadron. Therefore only one particle emerges on the hadron side of 



FIGURE 23. The non-zero contributions to the amplitude after color averaging. 



the collision. It may be possible to generalize these considerations to fully diffractive 
processes where the hadron is allowed to fragment. 
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